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ABSTRACT 



Context. The RAVE spectroscopic survey for galactic structure and evolution obtains 8400-8800 A spectra at 7500 resolving power at the UK 
Schmidt Telescope using the 6dF multi-fiber positioner. More than 300000 9</c<12 and |fo|>25° southern stars have been observed to date, 
kyl , Aims. This paper presents the first intrinsic examination of stellar spectra from the RAVE survey, aimed at evaluating their diagnostic potential 
j_j ■ for peculiar stars and at contributing to the general understanding of Luminous Blue Variables (LBVs). 

' Methods. We used the multi-epoch spectra for all seven LBVs observed, between 2005 and 2008, in the Large Magellanic Cloud (LMC) by 
the RAVE survey. 

Results. We demonstrate that RAVE spectra possess significant diagnostic potential when applied to peculiar stars and, in particular, LBVs. 
The behaviour of the radial velocities for both emission and absorption lines, and the spectral changes between outburst and quiescence states 
are described and found to agree with evidence gathered at more conventional wavelengths. The wind outflow signatures and their variability 
are investigated, with multi-components detected in S Doradus. Photoionisation modelling of the rich emission line spectrum of R 127 shows 
evidence of a massive detached ionised shell that was ejected during the 1982-2000 outburst. Surface inhomogeneities in the nuclear-processed 
material, brought to the surface by heavy mass loss, could have been observed in S Doradus, even if alternative explanations are possible. We 
also detect the transition from quiescence to outburst state in R 71. Finally, our spectrum of R 84 offers one of the clearest views of its cool 
companion. 



Key words. Stars: emission-line - Stars: winds, outflows - Stars: Wolf-Rayet - Magellanic Clouds - Surveys 



1. Introduction 

RAVE (RAdial Velocity Experiment) is an ongoing digital 
spectroscopic survey of stars in the magnitude interval 9< 
Ic <12, distributed over the whole southern sky at galactic lat- 
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itudes \b\ >25°. Spectra are recorded over the 8400-8800 A 
range, at a resolving power ~7500, with the UK Schmidt tele- 
scope feeding light to a spectrograph via the 6-degree Field 
(6dF) 150 fiber positioner. Via the determination of radial ve- 
locities, chemistries, temperatures, and gravities for a large 
number of high-latitude stars, the overarching science driver 
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for the survey is the investigation of the structure and evolution 
of the Milky Way. At the time of writing, RAVE has observed 
over 300,000 stars, Data Releases 1 and 2 have been published 
(Steinmetz et al. 2006; Zwitter et al. 2008), the third is coming 
(Siebert et al. 2009), and scientific exploitation of the data has 
begun (e.g. Smith et al. 2007; Siebert et al. 2008; Seabroke et 
al. 2008; Munari et al. 2008; Veltz et al. 2008). 

The essentially unbiased selection of targets guarantees that 
in addition to normal stars, peculiar ones are also observed. The 
type of peculiarities most easy to recognise with RAVE spec- 
tra are the presence of emission lines, absorption line splitting 
in double-line binaries, signatures of stellar winds, and com- 
bination spectra. The present paper is the first to explore the 
performance of RAVE in relation to the physics of peculiar 
stars, in particular, the Luminous Blue Variables (LBVs) of 
the Large Magellanic Cloud (LMC). As part of its tiling strat- 
egy, RAVE observed the region of sky occupied by the LMC 
in March 2005, March 2006, September 2006, and February 
2008. As LBVs are intrinsically very luminous, and those be- 
longing to the LMC are brighter than the /c<12 survey limit 
(indeed, many possess an HD identifier), they naturally enter 
the RAVE target set. A list of the LBVs observed by RAVE in 
the LMC is given in Table 1 . 

LBVs were first identified in M31 and M33 by Hubble & 
Sandage (1953) as extremely luminous hot stars that under- 
went irregular photometric variability of modest amplitude (a 
few mag) over timescales of years. Their bolometric magnitude 
generally exceeds -9.5, corresponding to a luminosity of the 
order ~10 6 L Q , close to their Eddington Limit. During quies- 
cence, the LBVs obey a tight temperature-luminosity relation, 
calibrated as log L/L G = 1.37 log T e g - 0.03 by van Genderen 
(2001). The existence of such a relation, their extreme intrin- 
sic brightness and easy identification (the latter provided by 
their variability), make LBVs important players in the cosmic 
distance ladder game (e.g. Davidson et al. 1989). The LBVs 
are sometimes also called Hubble-Sandage variables or S Dor 
variables from the archetype of the class, and are believed to 
be the immediate progenitor of Wolf-Rayet stars of the classic 
N-type (nitrogen rich and hydrogen poor). The LBV phase is 
not particularly long-lasting (~0.025 Myrs) and, together with 
the paucity of massive stars, accounts for the relative rarity of 
LBVs (Crowther 2007). The number of recognised LBVs dif- 
fers from author-to-author, ranging from 20 to 50 throughout 
the entire Local Group (e.g. van Genderen 2001; Weis 2003). 
The most famous LBVs include 77 Car and P Cyg, in the Galaxy, 
and S Dor in the LMC. 

After the main sequence phase, the most massive stars do 
not evolve all the way to the reddest part of the HR diagram 
but instead enter a phase of very high mass loss (10 4 — 10 -5 
M yr _1 ) and reverse their evolution back toward hotter effec- 
tive temperatures, becoming LBVs. The position of this turn- 
ing point depends upon the star's luminosity and sets what 
is called the Humphreys-Davidson Limit (cf. Humphreys & 
Davidson 1994). The LBVs lay to the hotter side of this T e g 
limit, which may be crossed during eruptive events of the most 
massive objects. The evolutionary sequence leading to the for- 
mation of LBVs depends on the mass of the progenitor. For 
stars more massive than 75 M it is O -> WN(H-rich) -» LBV 



Table 1. List of the LBVs observed by RAVE in the LMC, 
their date of observation and optical brightness at that time, and 
the radial velocities of the associated emission and absorption 
lines. 



HD 


name 


obs. 


date 


V 


RV em. 


RV abs. 
















(km /s) 


(km/s) 




35343 


S Dor 


2005 


03 


20 


9.55 


293.0 ± 1.4 


299.0 ±0 


.2 






2006 


03 


23 


9.45 


297.4 ± 1.0 


266.5 ±2 


.7 






2006 


09 


27 


9.60 


298.2 ± 1.9 


253.0 ±0 


.5 


269006 


R 71 


2006 


03 


23 


10.70 


199 ± 11 


186.7 ±2 


.1 






2006 


09 


27 


10.65 


193.4 ± 2.4 


205.4 ±2 


.2 






2008 


02 


04 


9.95 


197.5 ± 5.9 


157.7 ±2 


.6 


269445 


R 99 


2006 


03 


23 


11.50 


289.6 ± 5.1 










2006 


09 


27 


11.55 


288.4 ± 3.4 










2008 


02 


04 


11.60 


289.1 ± 4.0 






269929 


R 143 


2006 


09 


27 


11.15 




286.6 ±2 


.4 






2008 


02 


04 






249.1 ±2 


.6 


269216 


S 88 


2005 


03 


20 


10.40 




233.8 ±3 


.4 


269227 


R 84 


2006 


09 


27 


11.80 




254.2 ±1 


.2 


269858 


R 127 


2006 


09 


27 


10.40 


276.5 ±0.4 







-> WN(H-poor) -> WC -> SNIc, while for those of initial 
mass 40<M<75 M it is O ^ LBV WN(H-poor) -> WC 
— > SNIc (cf. Crowther 2007, and references therein). Direct 
spectroscopic evidence for LBV progenitors transitioning to 
supernovae is being obtained (eg. SN 2005gj - Trundle et al. 
2008). Only lower-mass O-type stars can experience a red su- 
pergiant phase just before or after the LBV phase (Szeifert et 
al. 1996; Smith et al. 2004), their evolutionary sequence being 
O -> LBV/RSG -» WN(H-poor) -> SNIb (where RSG=Red 
Supergiant). This is supported by the simultaneous presence of 
LBVs and RSGs observed in Westerlund 1, the most massive 
young open cluster of our Galaxy, with lxlO 5 M total mass 
and <25 M turn-off mass (Clark et al. 2005). To the best of our 
knowledge, this simultaneity is not known to occur elsewhere. 

The large amount of mass lost during the LBV phase is the 
critical stage that a very massive star must pass through before 
becoming a WR star. Mass ejection via major outbursts is far 
more efficient than that via steady winds, as most notably un- 
dergone by 77 Car during the 19th century when the star shed 
several M in less than a decade (Smith & Owocki 2006). The 
mass lost by the central LBV star frequently gives rise to cir- 
cumstellar nebulae (e.g. Langer et al. 1999, Weis 2003) that 
turn out to be enriched in processed material, such as nitrogen 
and helium, coming from the stellar interior. 

During the quiescent phase, the optical spectra of LBVs are 
characterised by effective temperatures ranging from 12 000 to 
30 000 K (B-spectral types) and strong emission lines of hydro- 
gen, Hel, Fell, Call and other singly ionised metals, often with 
P-Cygni profiles when observed at sufficiently high resolution 
(cf. Kenyon & Gallagher 1985). The mass loss rate in quies- 
cence is smaller and the wind faster. While not well-understood 
an LBV occasionally passes through an 'outburst' phase, dur- 
ing which the luminosity remains the same as in quiescence, 
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Fig. 1. Multi-epoch RAVE spectra of S Doradus. Paschen and NI (multiplets 1 and 8) absorption lines are indicated by comb 
markings. 
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Fig. 2. Multi-epoch RAVE spectra of R 99 and R 7 1 . 
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while in tandem the mass loss rate increases and the wind ve- the effective temperature is ~7500-8500 K regardless of the 

locity decreases. This pushes outward the pseudo-photosphere, LBV luminosity and temperature in quiescence, and the spec- 

with a consequent increase in the effective radius and decrease trum turns to that of A-type supergiants, while [Fell], Hel 

in the surface temperature (whose Wien peak moves from the emission lines, and P-Cygni absorptions weaken or disappear, 

ultraviolet toward the optical region). At outburst maximum, Wolf (1989; see Vink 2008 for an update) pointed out the ex- 
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istence of an amplitude-luminosity relation for LBV eruptions 
and called S Doradus instability strip the region of the HR dia- 
gram occupied by LBVs in quiescence. Stothers & Chin (1995) 
suggested the existence of a period-luminosity relation of the 
form M bo i = -12.5(±0.5) + 2.4(+0.5) logP where P, expressed 
in yrs, is the typical interval from one eruption to the next. 

2. RAVE Spectra of the LBVs 

The spectra that RAVE collected for the LBVs in the LMC are 
presented in Figs. 1-3, where the ordinate scale is kept con- 
stant for an easier inter-comparison. Differences in S/N among 
the multi-epoch spectra of the same star are mainly due to dif- 
ferences in the throughputs of different fibers and in the sky 
conditions. The radial velocities derived from the absorption 
and emission lines are given in Table 1 (for comparison, the 
heliocentric barycentric velocity of the LMC is +278 km/s, 
Mateo 1998, although significantly variable across the face of 
the galaxy). Table 1 also lists the V-magnitude at the time 
of the RAVE observation as estimated by consulting various 
databases of variable star amateur observers (AAVSO, VSNET, 
VSOLJ) and the ASAS Survey database (Pojmanski 2002). 

These observations indicate that the RAVE spectra are well 
suited to detect and study peculiar stars, and LBVs in particular, 
confirming early investigations in this wavelength interval by 
Munari (1998, 2002, 2003). The peculiarities stands out well 
in terms of absorption features, wind outflow signatures, and 
emission lines, and will now be considered in turn. 

2.1. Emission lines 

The ionised circumstellar gas is easily revealed by the rich 
emission line spectrum observable over the RAVE wavelength 
regime. The principal ions contributing emission lines at these 
wavelengths are HI, Hel, Call, NI, CII, CIH, OI, SI, [Fell], 
Fell, Fein, [CHI], and [C1III] . The emission lines identified and 
measured in the LBV RAVE spectra of Figs. 1-3 are listed in 
Table 2. They all belong either to OI, NI, [Fell], Hel, Call or to 
the Paschen series of hydrogen. 

The emission lines in RAVE spectra of LBVs are fairly 
"sharp". The expansion velocity derived from the width of the 
profiles (corrected for the FWHM(PSF)=40 km/sec instrumen- 
tal resolution) is on average 25 km/sec, which is the same 
as obtained from observations at more classical optical wave- 
lengths (Weis 2003). The only exception is R 99 (cf. Fig. 2), 
whose expansion velocity from Paschen lines in RAVE spec- 
tra is ~90 km/s, comparable to the Balmer line-based results of 
Notaet al. (1996). 

The intensity of the OI 8446 emission line can be boosted 
by fluorescence pumping from absorption of hydrogen Lyman- 
fi photons by OI in its ground state, as first pointed out by 
Bowen (1947). For the Lyman-/? fluorescence to be effective, 
the optical depth in Ha should be large, owing to the popula- 
tion of the n = 2 level by trapped Lyman-a photons. Such a 
fluorescence pumping is at work in R 127 (cf. Fig. 3), but not 
in R 99 (cf. Fig. 2), as derived by comparing the intensity of 
OI 8446 emission with that of the Paschen series of hydrogen, 



following Strittmatter et al. (1977) analytical modelling of the 
pumping mechanism. 

2.2. Wind outflows 

Significant mass loss from a central star is revealed in our LBV 
spectra by pronounced P-Cygni profiles. A selection of such 
profiles is presented in Fig. 4. The line showing a wind-driven 
P-Cygni profile most consistently is that of OI 8446 A. 

In a P-Cygni profile, the emission component traces the 
systemic velocity (eg. Castor & Lamers 1979), so that the blue- 
ward displacement of the absorption component relates to the 
velocity of the wind at the position in the expanding medium 
where the given line forms (Lamers & Cassinelli 1999). The 
RAVE spectra in Fig. 4 show that the wind velocity from OI 
lines is slow and quite similar in each of the surveyed LBVs, 
ranging from 130 to 210 km/sec. These velocities closely 
match those derived from high resolution optical and ultravi- 
olet spectra (eg. Leithereret al. 1992; Crowther & Willis 1994; 
Garcia-Lario et al. 1998). The RAVE spectra demonstrate that 
the velocity, intensity, and shape of the wind absorption com- 
ponent is highly variable for any given object, as illustrated in 
Fig. 4 by the OI profiles of S Dor and R 71 at three distinct 
epochs. For R 71, the wind absorption emerged only in the out- 
burst spectrum of 4 February 2008, contrary to the normal be- 
haviour of LBVs. S Dor showed significant variability of the 
OI absorption component during a period of constant bright- 
ness. The S Dor spectrum of 27 September 2006 shows two 
distinct wind components, at 140 and 210 km/sec. Wind ab- 
sorption profiles characterised by multi-components (variable 
with time) have been previously observed in high-resolution 
optical spectra of R 127 (Stahl et al. 1983) and AG Car (Stahl 
et al. 2001), extending over the same range of velocities de- 
tected in the RAVE spectra of S Dor. 

The high variability of the absorption components of P- 
Cygni profiles reflects the complexity of wind condition and 
structure in massive and LBV stars. The winds of LBVs are 
both aspherical and inhomogeneous, as revealed by the highly 
variable polarisation measured by Davies et al. (2005). The 
amount of mass loss (in M yr _1 ) in such stars relates to both 
luminosity (in L ) and effective temperatures (K) as log M = 
1 .738 log L - 1 .352 log r eff - 9.547, following the calibration of 
Lamers & Leitherer (1993). The wind velocity field is usually 
assumed to follow a yS-type velocity law v(r) = v 0O (l - R+Jff 
(cf. Schmutz et al. 1991). The smaller the f3, the shorter the 
distance the wind needs to travel before reaching the terminal 
velocity. Guo & Li (2007) found that 0.5<yS<0.7 fits the ob- 
servations of LBVs in quiescence, while p>\.5 is required for 
bright states. Thus, in quiescent LBVs, is reached within a 
very few stellar radii, beyond which the gas density declines 
as p(r) = p(r Q )[r /r] 2 , while during bright states v M is reached 
only at a much greater distance. 

A cursory inspection of the spectra of S Dor in Fig. 1 either 
(i) could suggest the presence of weak and variable emission 
cores within Paschen absorption lines, or (ii) could indicate 
that Call emission lines appear in emission within the adjacent 
Paschen 13, 15, and 16 lines, with no associated true P-Cygni 
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Fig. 3. RAVE spectra of R 127, S 88, R 84, and R 143. 
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absorption component. That this is not the case is illustrated by 
Fig. 5, where a "zoomed" view is provided for one of the S Dor 
spectra in Fig. 1. Figure 5 clearly illustrates how the the pres- 
ence of weak emission cores within Paschen absorption lines 
is actually due to the interplay with nearby NI absorption lines, 
and how Call lines do really present a P-Cygni absorption com- 
ponent not associated with the wings of adjacent Paschen lines. 

2.3. Absorption spectrum and chemical abundances 

The classification of the absorption spectrum of B-A-F stars at 
the RAVE wavelengths is straightforward, as illustrated by the 
observed MK spectral atlases of Munari & Tomasella (1999), 
Cenarro et al. (2001), le Borgne et al. (2003), Marrese et al. 
(2003), and Andrillat et al. (1995), as well as the LTE syn- 
thetic spectral atlas by Castelli & Munari (2001) and the atlas 
of non-LTE synthetic spectra of OB super-giants by Clark et 
al. (2005). The absorption spectra of LBVs can be classified 
using these atlases, with the caveat that both markedly non- 
solar surface abundance patterns and the pseudo-photosphere 
being placed in a moving, partially ionised wind, can confuse 
the picture (see Jaschek & Jaschek 1987 for a self-consistent 
definition of classification criteria attributes). Nevertheless, the 
distinction between quiescent (B-type absorption spectra) and 
outburst states (A-type absorption spectra), as well as the sig- 
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Fig. 4. Sample of P-Cygni profiles of LBVs from RAVE spec- 
tra. 
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natures of type la luminosity class are straightforward to recog- 
nise. 

Significant mass loss removes the external layers of mas- 
sive stars, exposing progressively deeper internal regions and 
reveals the mixing with nuclearly processed material, which 
manifests as a depletion of hydrogen, oxygen, and carbon, 
and an enhancement of helium and (especially) nitrogen (cf. 
Meynet & Maeder 2005). Hydrogen, helium, nitrogen, and 
oxygen possess strong lines within the RAVE wavelength in- 
terval, at the effective temperatures and gravities characteriz- 
ing LBVs. Actually, multiplets 1 and 8 over the RAVE range 
are the strongest NI features observable over the whole optical 
range (cf. Jaschek & Jaschek 1995). The hydrogen deficiency 
and nitrogen enrichment is quite obvious in the RAVE spectra 
of LBVs presented in this paper. There is no temperature, lumi- 
nosity, and metallicity combination, for which a normal stellar 
atmosphere with a solar abundance pattern can display NI lines 
stronger than Paschen ones as, for example, seen in the S Dor 
spectra of Fig. 1. This is both true for LTE and non-LTE syn- 
thetic spectra (Castelli & Munari 2001; Clark et al. 2005); for 
this to occur, it would require, at the same time, both an en- 
hancement in nitrogen and depletion in hydrogen. 

2.4. Radial velocities 

Table 2 presents the radial velocities of the LBVs measured 
from the RAVE spectra. The accuracy of the wavelength scale 
of RAVE spectra is -1.0 km/sec (Zwitter et al. 2008). The 
radial velocities of emission lines for the objects with multi- 
epoch observations are stable within a dispersion of 2.3 km/sec 
for S Dor, 2.4 km/s for R 7 1 , and 0.5 km/s for R 99. This is con- 
sistent with the long-term stability seen by Thackeray (1974) 
in LBV emission line radial velocities (within a few km/s), 
regardless of the phase (quiescent or outburst). Note that the 
observations of R 71 in Table 2 cover both the quiescent and 
outburst phases. 

The radial velocity of the emission lines in Table 2 are 
comparable to those derived from high-resolution optical spec- 
tra at various epochs. Wasselink (1956) derived +295 km/sec 
for S Dor, Thackeray (1974) found +195 km/s for R 71, and 
Weis reported +286 km/s for R 99. The situation for R 127 is 
less clear with a low accuracy value of +284 km/sec claimed 
by Feast et al. (1960), +292±9 km/s by Stahl et al. (1983), 
+267 km/s by Weis (2003), and +274+1 km/s by Stahl & 
Wolf (1986). The RAVE value of +276.5+0.4 km/s lies within 
2.5 km/sec of both the mean of these four determinations 
(+279+5 km/sec) and the most accurate value amongst them 
(+274+1), again consistent with a long-term constant value. 
Small differences in velocity of the emission component of a P 
Cygni profile are expected when the line strength is variable. 

On the contrary, the radial velocities of the absorption lines 
given in Table 2 are highly variable from epoch-to-epoch, by 
several tens of km/s. This relates to the variable positioning of 
the pseudo-photosphere within the accelerating portion of the 
out-flowing wind, driven at least in part by instability and in- 
homogeneities in the mass loss rate (Davies et al. 2005). This 
'breathing' pseudo-photosphere is also responsible for the low 



Table 2. Equivalent width of the emission lines identified and 
measured in the RAVE spectra of LBVs in the LMC. 



1 / i \ 

Ao(A) 


ion 


equivalent width (A) 






R 99 


R 127 


S Dor 






(27-09-06) 


(27-09-06) 


(20-03-05) 


8413 


P-19 


0.974 


0.556 




8438 


P-18 


1.167 


0.606 




8446 


[OI] 


0.216 


3.445 


0.148 


8451 


? 




0.825 




8467 


P-17 


1.270 


1.057 




8490 


[Fell] 


0.156 


1.023 




8498 
8502 


Call 
P-16 


}l.613 


0.460 


0.734 


8518 


Hel 




0.080 




8531 


Hel 


0.098 


0.107 




8542 
8545 


Call 
P-15 


j-2.065 


0.875 


1.138 


8583 


Hel 


0.269 


0.143 




8598 


P-14 


1.640 


0.725 




8615 


Hel 




0.326 




8617 


[Fell] 




0.068 




8629 


NI 


0.573 


0.188 




8636 


? 




0.101 




8662 


Call 


}2.851 




0.968 


8665 


P-13 


0.642 




8671 
8680 


Hel 
NI 




0.290 
0.392 




8683 


NI 




0.471 




8686 


NI 


} 1.377 


0.132 




8703 


NI 




0.240 




8712 


NI 




0.126 




8719 


NI 




0.090 




8729 


NI 




0.063 




8734 


Hel 


0.394 


0.122 




8750 


P-12 


1.232 


0.845 




8777 


Hel 


0.300 


0.177 








R 71 


S 88 


R 143 






(23-03-06) 


(20-03-05) 


(27-09-06) 


8446 


[OI] 


0.369 


0.242 


0.162 


8490 


[Fell] 




0.160 


0.142 


8617 


[Fell] 


0.743 




0.053 



amplitude, irregular photometric variability (with timescales 
of months to a few years) which is persistent in LBVs both 
in quiescence and outburst phases. The light-curves presented 
by Sharov (1975), Stahl et al. (1984), Humphreys & Davidson 
(1994), and Lamers (1995), amongst others, provide clear ex- 
amples of this restless variability. This variability can be further 
appreciated by comparing the RAVE values in Table 2 with the 
velocities given by Feast et al. (1960) for S Dor (+213 km/sec), 
R 71 (+198), R 143 (+263) and R 84 (+262). 

2.5. Photoionisation analysis 

The RAVE spectrum of R 127 provides a sufficient number of 
emission lines (cf Table 2) to attempt modelling of its photoion- 
isation structure. To this aim, we have used the CLOUDY code 
(Ferland et al. 1998) to evaluate the prevailing physical con- 
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Fig. 5. Expanded view of the S Doradus spectrum of Figure 1 
near the Paschen 12 + NI 8747, Paschen 13 + Call 8662, and 
Paschen 15 + Call 8542 features. The arrow points to the ab- 
sorption component of the Call P-Cygni profiles. 

ditions in R 127 at the time of the RAVE observations, and 
to assess the diagnostic potential of emission lines observable 
in the RAVE wavelength regime. Even though this exercise 
did provide consistent results (as outlined below), some de- 
gree of caution must be exercised, as the number of emission 
lines observable over the restricted wavelength range of RAVE 
is somewhat limited, and only a fraction of them is currently 
treated by CLOUDY (most notably, the NI and Hel lines are 
ignored). This type of modelling is best carried out by com- 
bining emission lines from the entire optical, ultraviolet, and 
infrared regimes, which provides access to a larger set of emis- 
sion lines from more ions and from a wider range of excitation 
conditions. Such an expanded set of input data could also allow 
exploration of more complex, multi-component emission envi- 
ronments, than the simple one-component geometry adopted 
here. 

RAVE optical fibers have an aperture of 7 arcsec, which 
corresponds to 1.7 pc at distance of the LMC. This is large 
enough to include a contribution (by unknown relative propor- 
tion) to the observed spectrum from the central star, its wind, 
the extended circumstellar nebula, any parent H II region, and 
any close optical companion (many LBVs are members of mas- 
sive and dense young stellar clusters). Concerning R 127, it is 
surrounded by a faint nebula 5 arcsec in diameter. Judging from 
spatially-resolved long slit spectra presented by Weis (2003), 
the light from this nebula (clearly visible on HST images ob- 
tained through an Ha+[NII] narrow-band filter) contributes a 
negligible fraction of the light compared with the unresolved 
central source. Similarly, there is is no parent HII region in 



which R 127 is embedded, so this is also not likely to be a prob- 
lem. The adaptive optic observations by Heydary-Malayeri et 
al. (2003) resolved R 127 into a dense cluster of stars, with 
two close optical companions to the LBV: one is 1 arcsec away 
and 5 mag fainter, and the other is 3 arcsec away and 2.7 
mag fainter. The spectrum of the latter star (as presented by 
Heydary-Malayeri et al. 2003) is deficient in emission lines, 
and consequently, is not contaminating the LBV emission line 
spectrum. In conclusion, only the ionised mass outflow from 
the LBV is contributing to the emission lines observed in the 
RAVE spectra. 

Walborn et al. (2008) have reviewed the recent photometric 
and spectroscopic history of R 127, and found it entered a long- 
lasting outburst phase in the early- 1980s, peaking in the V-band 
at 8.9 mag during 1988-1991, and then slowly declining toward 
quiescence (by the time of the RAVE observations in 2006, 
R 127 had reached reached this quiescent state). According to 
the analysis of Guo & Li (2007), the wind of LBVs in quies- 
cence is accelerated to its terminal velocity within a few stellar 
radii. Therefore, the vast majority of its mass is located exter- 
nal to the acceleration zone and is expanding at constant veloc- 
ity. In our modelling, we thus adopted p(r) = p(r Q )[r /r] 2 for 
the dependence of density on distance from the central star. To 
be conservative, we also explored other density dependencies 
(eg. wind accelerated to terminal velocity over significant dis- 
tances, density independent of distance from the central star), 
but abandoned them ultimately due to the resulting inferior re- 
sults. We assumed Case B nebular conditions (Osterbrock & 
Ferland 2006), spherical geometry with the photoionising star 
at the center, and chemical homogeneity throughout the gas. 
The stellar photosphere was assumed to emit as a black body. 
The chemical partition was taken to be solar for all elements 
except those allowed to vary. The initial value for the overall 
metallicity was taken to be [M/H]=-0.5. Having no observed 
emission line from carbon, and CLOUDY not treating the He 
and N lines recorded over the RAVE spectral range, we as- 
sumed (and kept fixed) their abundance to the average value 
found in other LBVs - i.e., a lOx enhancement for N and 2.5x 
for He, and a 0.2x depletion for C (cf. Crowther & Willis 1994, 
and references therein). These fixed values for He, N, and C are 
close to those found by Lamers et al. (2001) for the external, 
spatially resolved, nebula around R 127. 

The emission lines that we used in the modelling, along 
with their observed and computed flux ratios, are shown in 
Table 3. The difference between computed and observed fluxes 
is ~10%. This can be regarded as an excellent match, reflecting 
the limitations/capabilities of such photoionisation modelling, 
regardless of the object under analysis (eg. novae, planetary 
nebulae, galactic nuclei, etc.) 

The output parameters of the modelling are given in 
Table 4, where £ is the filling factor (the ratio of filled-to-total 
volume in the ejecta). The derived temperature, radius, and lu- 
minosity are consistent with other modelling attempts in the 
literature, despite being based on completely different sets of 
input data. For example, Stahl et al. (1983) derived from optical 
and IUE high resolution observations of R 127 (during the qui- 
escence phase of the early-1980s), r eff =16000 K, R=150 R , 
and Mb o i=-10.6. The hydrogen and oxygen under-abundances 
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Table 3. Comparison between computed and observed emis- 
sion lines. Intensities are scaled to Paschen 14. 



Ao(A) 


ion 


relative flux 






observed 


computed 


8446 


[01] 


4.75 


4.73 


8498 


Call 


~0 


<0.01 


8542 


Call 


~0 


<0.01 


8545 


P-15 


1.21 


1.61 


8598 


P-14 


1.00 


1.00 


8617 


[Fell] 


0.09 


0.08 


8662 


Call 


~0 


<0.01 


8665 


P-13 


0.89 


1.14 


8750 


P-12 


1.17 


1.29 



we have obtained are similar to those found in other LBVs 
(Lennon et al. 1994; Crowther & Willis 1994; Venn 1997; 
Smith 1997), in R 127 itself (Lamers et al. 2001), and expected 
by theoretical evolution models of massive stars that include 
rotation (Meynet & Maeder 2005). 

The results in Table 4 indicate that the emission lines of 
R 127 originate in an ionised shell detached from the central 
star, extending from 19 to 27 stellar radii. The mass within this 
ionised shell is: 

M she u = I f ~ \nr 2 p H {r)dr = 1.33 x 10~ 3 M (1) 

where X is the hydrogen mass fraction, and p H (r) = p^(r m /r) 2 . 

The existence and meaning of such a detached shell de- 
serves comment, as fits squarely into the coherent picture 
emerging as a result of the full suite of data associated with 
R 127. 

Weis (2003) and Smith et al. (1998) determined from high 
resolution spectra that the bulk velocity of the mass outflow 
during the recent outburst was ~10 km/sec. Material expand- 
ing with this velocity should have left R 127 during 2000, in 
order to have reached (by the time of our RAVE observations) 
a distance from the central star equal to that of the inner border 
of the ionised shell in our model. It is interesting to note that the 
lightcurve of Walborn et al. (2008) shows R 127 started a rapid 
decline in brightness in 1999 that ended in 2000 when the star 
reached its "stable" quiescent brightness. Thus, we argue that 
the end of the outburst in 2000 brought with it also the cessa- 
tion of an enhanced mass loss phase, producing a cavity around 
the star that was only tenuously filled by the much lighter wind 
characterizing the subsequent quiescence conditions. 

The material at the external ionisation boundary of the shell 
(17.6 AU), expanding at 10 km/sec, should have left R 127 
2.5 yr before the material at the inner edge. The mass loss rate 
necessary to produce the observed shell is therefore 

M = ^ = 5 x 10- 4 M /yr (2) 

This is a magnitude larger than the mass loss rate derived by 
Stahl et al. (1983) for R 127 before the outburst. A large in- 
crease in the mass loss is a characteristic feature of LBVs dur- 
ing outburst phases (Humphreys & Davidson 1994). The ma- 



Table 4. Basic parameters for the central star and ionised cir- 
cumstellar material derived by photoionisation modelling of 
the emission line spectrum of R 127. 
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Fe/Fe 
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Ca/Ca 
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N/N 


(10) 
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terial lost earlier in the outburst had, by the time of the RAVE 
observations, moved exterior to the ionisation boundary, where 
the irradiation from the central star is no longer effective in ion- 
ising the gas, and therefore it was no longer contributing to the 
emission line spectrum. 

This scenario is well matched by the spectroscopic evolu- 
tion that Walborn et al. (2008) presents for R 127. The 1999- 
2000 timeframe, during which the star ejected the material 
comprising the ionised shell seen in the RAVE spectra, was 
characterised by the disappearance of P-Cygni absorption pro- 
files and by the appearance of very sharp emission line profiles. 
These are indicative an extended atmosphere moving outward 
at low velocity. Later spectra, when the significant mass loss 
had finished and the detaching of the massive shell had ini- 
tiated, started showing a progressive broadening of the emis- 
sion lines and re-emergence of P-Cygni profiles characteristic 
of a less-significant and faster mass loss regime. The winds 
of LBVs are line-driven, and the drop in outward velocity and 
large increase in mass loss rate during outburst phases is caused 
by the driving effect of Felll lines below the sonic point as the 
effective temperature decreases (Vink et al. 1999). The bifur- 
cated wind regimes between quiescence and outburst phases, 
referred to as the 'bi-stability' jump (Pauldrach & Puis 1990), 
is a defining characteristic of the S Dor instability strip (Vink 
& de Koter 2002). 

3. Remarks on Individual Objects 

3.1. S Doradus 

The three RAVE spectra presented in Fig. 1 were obtained 
while S Dor was stable during its outburst, at V~9.5 mag. Much 
like the Thackeray (1965) analysis, the continuum of the RAVE 
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spectra of S Dor can be readily classified as A2/3 Ia + e, by com- 
parison with the Munari & Tomasella (1999) atlas. 

Quite interesting is the reversal in the ratio of equivalent 
width (REW) of hydrogen and NI absorption lines that oc- 
curred between 20 March 2005 and 27 September 2006 spectra 
(cf. Fig. 1). The ratio changed from REW=0.7 to REW=3.0. 
This happened at stable brightness, surface temperature, and 
wind outflow velocity, while the intensity of the Call emission 
components increased by 1.6x between the two dates. It would 
be tempting to interpret the reversal of REW as an indication of 
chemical inhomogeneities in the nuclearly processed material 
progressively exposed at the stellar surface by the on-going sig- 
nificant mass loss. To support this interpretation additional and 
higher resolution observations are however necessary. In fact, 
the observed variation in the intensity of the Call emission lines 
could suggest the presence of unresolved and variable emission 
cores as a (contributing) reason for the observed changes in the 
equivalent widths of nitrogen and hydrogen lines. 

3.2. R 71 

The LBV nature of R 71 was discovered by Thackeray (1974), 
who reported that the star oscillates between 9.8 and 10.9 mag 
(see also Lamers 1995) and that the corresponding spectral 
types ranged from Aleq to B2.5Iaep. During bright optical 
phases, Balmer lines show deep P-Cygni profiles and the [Fell] 
emission lines disappear. 

The first two RAVE observations of R 71 were obtained 
when the star was in quiescence at V~10.7, while the third 
was during outburst, at V=9.95. The corresponding spectral 
changes in Fig. 2 bear a striking resemblance to the behavior 
described by Thackeray (1974). The RAVE quiescence spectra 
are clearly those of an early B super-giant (strong Hel absorp- 
tion lines, no OI, and broader Paschen lines) with superimposed 
strong [Fell] emission lines. The outburst spectrum is that of an 
early A super-giant, where the Hel lines have disappeared and 
Call and OI have appeared, along with the reinforcement of 
NI and the sharpening of Paschen lines. Similar to the optical 
spectra, in the outburst spectrum recorded by RAVE, the [Fell] 
has nearly disappeared and OI has developed a strong P-Cygni 
absorption component. 

3.3. R 84 

R 84 is an enigmatic object, suggested to be a dormant LBV 
by Crowfher et al. (1995). Walborn (1977) classified the blue 
optical spectrum as Olafpe, with many strong emission lines 
(including Hell 4686) and P-Cygni profiles. Higher resolution 
blue optical spectroscopy by Stahl et al. (1985) show the same 
characteristics. Breysacher et al. (1999, and references therein) 
summarised the determinations of the physical parameters ob- 
tained by various authors, which cluster around a temperature 
of 30000 K, a radius of 30 R , a luminosity 6xl0 5 L Q and a 
mass loss rate 3xl0~ 4 M /yr. 

The spectral appearance at longer wavelengths is however 
in conflict with that at the bluest end of the optical range. 
TiO absorption bands were observed shortward of 6000 A by 



Sanduleak & Philip (1977) on objective prism plates, and by 
Cowley & Hutchings (1978) in slit spectra. Allen & Glass 
(1976) found no evidence of them in far red spectra, but re- 
ported that the Mgl-b triplet could be present weakly in ab- 
sorption and that the IR colors implied the presence of a cool 
super-giant. Stahl et al. (1984) noted that the IR colour excess 
was too large to result only from a late-type companion, and 
suggested the presence of dust in R 84, which was confirmed 
by the mid-IR observations of Glass (1984). 

These apparently contradicting reports on the spectral ap- 
pearance support a binary or optical pair nature of R 84. Some 
photometric variability of one or both component stars is the 
probable cause of their brightness ratio to change with time. 
Schmutz et al. (1991) argued strongly in favor of the LBV and 
the cool star being merely along the same line of sight and not 
members of a binary system. It should be noted that adaptive 
optic observations by Heydari-Malayeri et al. (2003), which 
have a resolution limit of 0.12 arcsec (corresponding to 6000 
AU at the distance of the LMC), did not resolve the central 
star. Stahl et al. (1984) found a good fit to the IUE-optical-IR 
spectral energy distribution by combining an M4Ia spectrum 
with a BOIa spectrum, and with the addition of a circumstellar 
dust component. The spectral type of the cool super-giant was 
given as ~M2 by Cowley & Hutchings (1978). 

The RAVE spectrum of R 84 in Fig. 3 presents the clearest 
view so far published of the cool super-giant. It has no emis- 
sion lines superimposed and a x 2 fitting to the same Munari 
et al. (2005) synthetic spectral library used in the automatic 
analysis of RAVE spectra suggests atmospheric parameters of 
r eff =3925+155, [M/H]=-0.39±0.26, logg=2.87±0.91, and a 
heliocentric radial velocity of 254.2+ 1 .2 km/sec. The radial ve- 
locity reported by Schnurr et al. (2008) is 255 km/sec. The rela- 
tively large errors of the x 2 fitting are due to the contamination 
from the superimposed hot super-giant spectrum, whose main 
effect is to provide an overall veiling of the absorption lines. 
This has a particularly adverse effect on the wings of the Call 
lines (most sensitive to log g), amplified by the perturbation of 
the shape of the adjacent continuum by the wide underlying 
Paschen absorption lines of the hot component. The^ 2 temper- 
ature would suggest a spectral type around K7, with the metal- 
licity being appropriate for massive stars in the LMC. 

4. Conclusions 

The primary science driver for the RAVE survey is the in- 
vestigation of the structure and evolution of the Milky Way. 
Nevertheless, as a natural byproduct of the survey, the reso- 
lution and diagnostic potential of the associated stellar spectra 
ensures a valuable resource for the stellar astrophysics commu- 
nity. 

The present paper is the first in our series to explore the 
performance of RAVE in relation to the physics of peculiar 
stars, in particular, the Luminous Blue Variables of the Large 
Magellanic Cloud. The RAVE spectra provide a clear and com- 
prehensive view of the LBVs, highlighting their temporal evo- 
lution and star-by-star differences. In particular, these have: (i) 
extended the temporal baseline and wavelength domain over 
which the radial velocities of emission lines have been ob- 



10 



U.Munari et al.: RAVE spectroscopy of LBVs in LMC 



served to remain constant, which reinforces the notion of an 
origin in a spherically symmetric component of the circumstel- 
lar environment, dynamically de-coupled from the central star; 
(ii) mapped the large variability of the absorption lines (both 
in radial velocity and profiles) in response to the "breathing" 
of the pseudo-photosphere forming in the heavy, out-flowing 
wind; {Hi) documented the high variability and multi-shell na- 
ture of the wind P-Cyg absorption profiles; (z'v) provided a clear 
recipe of the spectroscopic changes characterising the transi- 
tion of LBVs from quiescence to outburst phases and vice- 
versa, in a close match to evidence gathered at more conven- 
tional (blue) optical wavelengths; (v) offered the opportunity to 
test quantitatively the soundness of photoionisation modelling 
of the rich emission line spectrum observable over the RAVE 
range; (vii) discovered in R 127 the presence, and quantified the 
physical properties of, a massive detached ionised shell which 
was ejected during the 1982-2000 outburst. 
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